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Summary 

Flame temperatures in nonsooting and sooting environments 
were successfully measured by radiometry for premixed 
propane-oxygen laminar flames stabilized on a water-cooled, 
porous sintered-bronze burner. The measured temperatures 
in the nonsooting flames were compared with fine-wire 
thermocouple measurements. The results show excellent 
agreement below 1700 K, and when the thermocouple 
measurements were corrected for radiation effects, the 
agreement was good for even higher temperatures. The 
benefits of radiometry are (1) the flow is not disturbed by an 
intruding probe, (2) calibration is easily done by using a 
blackbody source, and (3) measurements can be made even 
with soot present. The theory involved in the radiometry 
measurements is discussed, as well as the energy balance 
calculations used to correct the thermocouple temperature 
measurements. 


Introduction 

Combustion gas temperatures are difficult to measure since 
radiation losses caused by the high temperatures associated 
with flames produce errors in the thermocouple readings. 
Temperature measurements are also limited by the ther- 
mocouple material. In addition, a special problem arises in 
measuring combustion gas temperatures when soot is present. 
The soot quickly builds up on the thermocouple, continuously 
changing the emissivity of the thermocouple, which makes 
radiation heat transfer corrections very difficult. 

Engine designers are interested in soot research since the 
presence of soot increases the radiation heat transfer to engine 
walls, reducing engine life. Interest in soot formation in flames 
led to a desire to nonintrusively determine the soot volume 
fractions, particle sizes, and number densities using laser 
scattering-extinction in both diffusion (refs. 1 and 2) and 
premixed (refs. 3 and 4) flames. Local gas temperature also 
is an important factor in the production of soot. To further 
understand the soot formation and burnup processes, it is 
necessary to measure the temperatures in flames and in the 
combustion gases downstream of the flame fronts. A viable 
measurement technique is infrared radiometry. This method 
is nonintrusive and does not require seeding the flow with any 
foreign particles or chemicals. The infrared monochromatic 


radiation method of gas pyrometry was first developed by 
Schmidt (ref. 5) in 1909. Much later, this technique was used 
as the basis for the modern, commercial infrared radiometer 
which was used in this study. 

This investigation determined temperature profiles in both 
nonsooting and sooting premixed propane-oxygen laminar 
flames. Infrared radiometry was used, and the results in the 
nonsooting flames were compared to fine-wire thermocouple 
measurements. The premixed flames were stabilized on a 
water-cooled, porous sintered-bronze burner (fig. 1). This 
burner produced a very uniform flow field in which the 
equivalence ratio and cold gas velocity could be varied with 
precise micrometer valves. This type of burner has been used 
in many research programs (refs. 6 to 10). 

Theory 

The laws of radiation originated by Stefan-Boltzmann, Wien, 
Planck, and others establish definite and precise relationships 
between an object’s temperature and composition and the 
radiation it emits. Conversely, by making accurate spectral 
measurements of the radiation received from a remote object, 
much can be learned about its thermal history and composition 
(ref. 11). These spectral measurements are made by using 
radiometry. The conventional process of radiometric data 
interpretation uses the basic laws of radiation together with 
the detailed characteristics of the optical, electronic, and 
detector elements of the instrument. 

One approach treats the radiometer as a system which has 
a transfer function, or response factor, K s , which converts the 
radiometer output voltage and scale settings to an absolute 
radiation measurement. The radiometer system response factor 
is defined as (ref. 12) 


where 

S instrument electronic attenuation factor 
V recorded output voltage 
N T effective radiance, W/cm 2 -sr 

N r effective radiance of radiometer internal reference, 
from calibration, W/cm 2 -sr 
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Figure 1.— Cross section of flat flame burner. 


(A complete list of symbols is provided in appendix A for the 
reader’s convenience.) First, A r must be determined by using 
a graphical technique along with the calibration data for the 
radiometer; N T is defined as 

N t = \ Xl R x N x d\ (2) 

% 

where 

R x normalized product of detector response and optical 
filter transmission characteristics (provided by 
instrument calibration) 

N x blackbody spectral radiance at selected temperature 
given by Planck’s equation: 



where 

Cj first radiation constant, 1.19X10 4 W-pm/sr-cm 2 
C 2 second radiation constant, 1.438 xlO 4 pm-K 
T blackbody temperature, K 

X all wavelengths in range of detector and filter (from 
Xj to X 2 ) 

In order to determine the radiometer detector must be 
calibrated over a range of wavelengths, and the detector 
response must be normalized and plotted versus wavelength 
(fig. 2). Then the optical filter must be calibrated in a similar 


fashion to produce a normalized calibration curve as shown 
in figure 3. The values of the detector response and the filter 
response at a given wavelength are multiplied together, and 
the resulting curve is again normalized. This new curve is R x , 
which is used in equation (2) along with equation (3) to 
determine N T . 

Evaluating the integral in equation (2) for a selected 
temperature will give its corresponding value of N r . These 
results were tabulated so that in the future, when an experimen- 
tal value of N t is found, the table will provide the 
corresponding temperature (ref 12). For instance, the table 
is immediately useful to determine N r for equation (1). The 
internal reference temperature was measured at the time of 
calibration. We now use the radiance versus temperature (N T 
versus T) table to look up this reference temperature and find 

A In 
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Figure 3.— Radiometer filter response. 
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After using a blackbody source to supply a known tempera- 
ture to calibrate our system and determine K s , we can proceed 
to determine unknown temperatures by using equation (1) 
rearranged as 

SV 

N T = N r + - (4) 

K s 

Experimental Method 

When attempting to measure a flame temperature, we cannot 
assume that the flame is a blackbody (emissivity = 1). The 
flame emissivity must be determined experimentally. 
Neglecting the effects of the surrounding atmosphere and 
assuming the flame is at a uniform temperature, we seek the 
flame spectral emittance e X i/ at \ 1? which is chosen to be 
4.4 pm for our study in order to collect C0 2 radiation, since 
C0 2 is plentiful in combustion products. The following 
identity will be used: 


V f 2 voltage with flame plus transmitted blackbody 
radiation 

Two equations can be written for cases 3 and 4 seen in 
figure 4: 

Vf\ ~ K r \\f+ ot\x i/V\ibf ( 5 ) 

and 

Vfl = V bb r Mf + <*\lfV\\bf ( 6 ) 

where V X \ b f is the voltage that would be recorded if the flame 
were a blackbody. Since the flame is not a blackbody, this 
voltage (radiation intensity) must be multiplied by the flame 
absorptance in these equations. 

Subtraction of equations (5) and (6) yields 

Vf2 “ Vf\ — T X\ f(Vbb “ K) 


a xi / - e xi f 


which can be solved for the flame transmittance, 


where 

a xi/ spectral absorptance of the flame 
e xi / spectral emittance of the flame 

and we can define r X i/to be the spectral transmittance, such 
that 1 = a + r. 

Since the voltage measured by the radiometer is a measure 
of the radiation intensity, voltages and radiation intensity will 
be considered to be interchangeable in the following analysis. 
The method (ref. 5) requires measurement of four voltages 
(fig. 4): 

V 0 voltage with blackbody covered with shutter 
V bb voltage with blackbody radiation only 
V fl voltage with flame radiation only 


T W = 



( 7 ) 


Now, either equation (5) or (6) can be used to find V Mb f, 
the voltage that would be measured if the flame were a 
blackbody source radiating at the flame temperature. For 
example, using equation (5) yields 


Vx W = 


Vf i ~ Krxif 

a Xi f 


which, upon using a = 1 — r, becomes 
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Figure 4.— Required voltage readings for radiometer flame temperature 
measurement. 


./ - V f ! “ V o r Wf /o\ 

v \\bf : W 

1 - T\l f 

By returning to the radiometer equation (4) and using the 
voltage V \lbf which corrects the flame to a “black” flame, 
a value for N T is determined which is the effective flame 
radiance corrected for flame emittance. 

A f =n (9) 

K s 

By using this N T , the corresponding flame temperature is 
found from the table. A curve fit to this table is used in a 
computer program for calculating a temperature from the four 
measured voltages. This program is shown in appendix B. It 
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is written in the APL language and can be run on IBM-PC/AT 
or compatible using the APL* Plus System (ref. 13). 

With this method, if there is soot in the flame, the soot will 
radiate in the 4.4-ptm wavelength along with the C0 2 . It is 
unnecessary to differentiate between the radiation from the soot 
and the radiation from the C0 2 , however, because the C0 2 
and soot are approximately at the flame temperature. The 
radiometer receives the radiation from both and converts it 
to the flame temperature. Of course, the flame emissivity will 
be different when there is soot, but the emissivity is measured 
and used in the calculations when determining the flame 
temperature. 

Thermocouple Temperature 
Measurements 

An energy balance between the convection to and the net 
radiation away from the spherical thermocouple bead is 


required to extract the true gas temperature from the bead 
temperature reading. Because the leads are aligned along an 
isotherm, we assume that there is no conduction down the 
leads. Figure 5(a) shows a photograph of the thermocouple 
apparatus. The upper wire is platinum/platinum-rhodium 
for flame temperature measurements, and the lower wire 
(Chromel/Alumel) can more accurately determine the lower 
temperatures of the burner surface (near 800 K). Figure 5(b) 
shows the position of the thermocouple bead in the flame. 

The convective heat transfer coefficient is approximated 
by using Nusselt number Nu correlations, where for a sphere 
in the creeping flow regime (i.e., Reynolds number Re 
approaching zero) (ref. 14) 

Nu = — = 2 + 0.37 Re 06 Pr 033 (10) 

k 

for all Reynolds and Prandtl Pr numbers of interest in the flat 
flame flow field. Here h is the average heat transfer coefficient 




Figure 5 . —Thermocouple used in flame temperature measurements. 
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(W/m 2 -K), D is the bead diameter (m), and k is the gas 
thermal conductivity (W/m-K). For the propane-oxygen flame 
(assuming only water and carbon dioxide products), for 1 mole 
of propane, 4 moles of C0 2 and 3 moles of H 2 0 are formed. 
The thermal conductivities of C0 2 and H 2 0 are formed as a 
function of temperature (ref. 15), and then the thermal 
conductivity for the mixture can be calculated as 


over a circular plate (with a radius of 3 cm), giving 
F bp =F bB = 0.276; then F boo = 1 — F bp - F bB = 0.448, 
while the gas surrounding the bead gives F bg 1. There was 
no soot present in the flame for the thermocouple meas- 
urements. This reduces equation (11) to 


<lr =W\ 


7/ - 0.08287/ - 0.02767/ - 0.03457/ - 0.4487/ 


k = 



3 , 

+ - k H2 o 


This radiative heat transfer is balanced by the convective heat 
transfer q ” c = h(T g - T b ) - hAT — q r \ so the gas temperature 
is found from 


An energy balance is then made around the thermocouple 
bead by using an optically thin approximation, where the gas 
and soot particles emit but do not attenuate radiation. Radiation 
is exchanged between the bead and the burner surface, the 
upper stabilizing plate, the hot gases, the soot, and the ambient 
atmosphere. Surfaces are assumed to be in local thermal 
equilibrium, diffuse and gray, however, the ambient atmos- 
phere is assumed to be black. The net radiation flux q" out 
of the bead is then 

q" r = e b oT A - a b F bp e p oT p - a b F bB e B < jT$ - a b F bs e s oT ’/ 

- a b F bg e g oTg - a h F bx oTi (11) 

where the subscripts indicate— fc, bead; /?, upper plate; B , 
burner surface; 5, soot; g , gas; and oo, ambient; and where 
e is the emissivity, a is the Stefan-Boltzmann constant, a is 
the absorptivity, and Fy are view factors. 

Kirchhoffs law and the gray assumption are used to 
substitute the emissivity e b for the absorptivity a b of the bead. 
This emissivity is found as a function of temperature from 
(ref. 16). 

e b = 1.507X 10 _4 (7 A ) — 1.596X10 _8 (7* 2 ) (12) 

for 0 < T b < 2230 K. 

The emissivities of the steel plate and bronze burner sur- 
face are approximately (ref. 17) e p * 0.3 and e B ~ 0.1. The 
emissivity of the gases are found from 

e g ~ e co 2 + e H 2 o ~ Ae (13) 

where the e of each of the gases is dependent on the partial 
pressure and the path length through the gas from the bead. 
Figures showing the emissivities for carbon dioxide and for 
water vapor are found in Siegel and Howell (ref. 17). The 
Ae, which corrects for the spectral overlap of the absorp- 
tion bands of C0 2 and H 2 0, was found to be negligible 
for this small flame (a path length of 15 mm). Then, 
e* = 0.025 + 0.0095 = 0.0345. 

The shape factors were determined by using a model (ref. 
17) of a sphere (the thermocouple bead was nearly spherical) 


Tg = T b + e -^- [7/ - 0.08287/ - 0.02767/ 

- 0.03457/ - 0.4487/] (14) 

Iteration is required since the gas conductivity is a function 
of temperature and this equation includes a fourth power term 
in T g . For example, for a thermocouple reading of 1780 K, 
the burner surface was nearly 800 K, the plate was approxi- 
mately 1000 K, and the gas temperature was 

T g = 1780 + 58 - 0.47 - 0.07 - 2.42 - 0.02 = 1835 K 

This shows that the major contributions to the thermocouple 
correction come from the radiation loss from the bead and from 
the additional radiation from the hot flame gases. The 
contributions from the burner surface, the plate, and the 
ambient atmosphere were negligible. An error analysis of 
this method (ref. 14) showed that, in a nonsooting flame, AT 
(=T g — T b ) was ± 15 percent accurate with the errors in T b 
and Tg both approximately ±3 percent. 

Results 

The flame temperature measured 15 mm above the burner 
surface and determined by the radiometer is shown as a 
function of equivalence ratio in figure 6, for three cold gas 
velocities— 6, 9, and 12 cm/s. These flames were nonsooting 
so that they could be compared with fine-wire (0.08-mm- 
diameter) platinum/platinum-rhodium thermocouples. For the 
higher flow rates (9 and 12 cm/s), there is an increase in 
temperature with an increase in equivalence ratio until the ratio 
is nearly stoichiometric. The temperature dips near the 
stoichiometric equivalence ratio, then continues to climb until 
near an equivalence ratio of 1 .3, where the temperature begins 
to drop again. This is not the trend for adiabatic well-stirred 
reactor combustion of propane-oxygen mixtures as a function 
of equivalence ratio. 

Figure 7(a) shows the expected temperature profile predicted 
by equilibrium chemical kinetics calculations (ref. 18). Notice 
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Figure 6.— Flame temperature determined by radiometry. 

the peak in temperature near the stoichiometric equivalence 
ratio. Figure 7(b) is a plot of the laminar flame velocities 
as a function of equivalence ratio for propane-air flames 
(ref. 19). The peak velocity is also seen to occur near the 
stoichiometric equivalence ratio. The velocity will be much 
higher when pure oxygen is used instead of air, as shown in 
the relative plot of maximum flame velocities as a function 
of mole fraction of oxygen shown in figure 7(c) (ref. 20). So, 
both the velocity and temperature profiles are expected to peak 
near the stoichiometric equivalence ratio. This can be used 
to explain the dip in flame temperature near the stoichiometric 
equivalence ratio for the flat flame burner. The flame speed 
increases with increasing equivalence ratio, and the flame 
begins to nestle closer to the water-cooled (i.e., nonadiabatic) 
burner surface, losing heat to the burner. Since the flame speed 
is maximum at the stoichiometric equivalence ratio, the 
greatest heat loss to the burner is seen here, even though 
the flame temperature should be the highest. Above the 
stoichiometric equivalence ratio, the flame speed begins to 
decrease, allowing the flame to lift farther from the burner 
surface, causing it to lose less heat to the water-cooled burner. 
This lessening of the heat transfer to the burner overcomes 
the lowering of the flame temperature due to the equivalence 
ratio. This effect is seen as an overall increase in the flame 
temperature until the equivalence ratio effect predominates 
when the flame has reached a position far enough away from 
the burner to be essentially adiabatic. At this point (near an 
equivalence ratio of 1 .3), the temperature begins to drop with 
an increase in equivalence ratio as expected. 

The low flow rate case (6 cm/s) has a flow velocity that 
allows the flame to remain at nearly the same position for all 
equivalence ratios. The temperature increases with increasing 
equivalence ratio to near the stoichiometric equivalence ratio 
and decreases thereafter as predicted by well-stirred reactor 
calculations. 



(a) Equilibrium flame temperatures for propane-oxygen premixed flames. 

(b) Laminar burning velocity as function of equivalence ratio for propane-air 
flames. 

(c) Effect of oxygen concentration on relative flame velocities using propane 
fuel. 

Figure 7.— Predicted flame temperatures and laminar burning velocities for 
adiabatic, premixed flames. Initial mixture temperature, 298 K; all flames 
at 1 atm. 

Figure 8 shows a comparison of fine-wire thermocouple 
readings with the radiometer temperature measurements of 
figure 6. The 45° line shows where perfect agreement would 
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Figure 8.— Comparison of radiometer and thermocouple temperature 
measurements. 


Figure 9.— Soot volume fractions and flame temperatures as functions of 
equivalence ratio. Height above burner, 15 mm; cold gas velocity, 3.5 cm/s. 


occur. The two techniques agree very well for temperatures 
below 1700 K. Above this temperature, radiation losses from 
the thermocouple begin to have a greater effect and the 
thermocouple reads a lower temperature. The thermocouple 
readings for the 12 cm/s case were corrected, and the corrected 
results are seen as solid symbols in figure 8. 

The radiometer was then used to measure temperature in a 
sooting flame (0 > 1.7) as a function of equivalence ratio 0, 
for a cold gas velocity of 3.4 cm/s, 15 mm above the burner. 
This temperature is shown in figure 9 along with the soot 
volume fraction (fraction of flame volume filled by the soot 
particles in parts per million) measured by a scattering- 
extinction laser diagnostic technique (ref. 4) for the same 
flame. As expected, the flame temperature decreases with 
increasing equivalence ratio, for equivalence ratios over 1.0, 
partially due to the increase in soot particles radiating heat 
to the environment. 


Conclusions 

Flame temperatures in both nonsooting and sooting environ- 
ments have been successfully measured by radiometry. In the 
nonsooting flames, comparisons with fine- wire thermocouples 
showed excellent agreement below 1700 K, and when the 
thermocouples were corrected for radiation effects, the 
agreement was good for even higher temperatures. The 
benefits of radiometry are (1) the flow is not disturbed by an 
intruding probe, (2) calibration is easily done by using a 
blackbody source, and (3) measurements can be made even 
with soot present. 


Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, November 28, 1988 
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Appendix A 
Symbols 


Q,C 2 first and second radiation constants, respectively 
D bead diameter 

Fy view factors 

h average heat transfer coefficient 

K s transfer function (or response factor) 
k gas thermal conductivity 

N r effective radiance of radiometer internal reference 
(from calibration) 

N t total effective radiance 

N x blackbody spectral radiance at selected 

temperature 

Nu Nusselt number 

Pr Prandtl number 

q” convective heat transfer 

q" radiative heat transfer 

R x normalized product of detector response and 

optical filter transmission 

Re Reynolds number 

S instrument electronic attenuation factor 

T blackbody temperature 

V recorded output voltage 

V\\ b f voltage that would be recorded if flame were a 
blackbody 

a absorptivity 

a \if spectral absorptance of flame 


e emissivity 

e\i f spectral emittance of flame 

X all wavelengths in range of detector and filter 

X lf X 2 minimum and maximum detectable wavelengths, 
respectively 

a Stefan-Boltzmann constant 

r X i/ spectral transmittance such that 1 = a 4- t. 

<t> equivalence ratio 

Subscripts: 

B burner surface 

b bead 

bb blackbody 

c convective 

/ 1 flame radiation only 

/2 flame radiation plus transmitted blackbody 

radiation 

g gas 

o blackbody covered with shutter 

P upper plate 

r radiative 

5 soot 

X wavelength 

oo ambient 
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Appendix B 

Program for Calculating Flame Temperature 


I 


V X C02TEMP Y ; ALN ; AT ; EMM ; NTBB ; NTF ; VOLT ; NREF ; SAT ; TAU 
[1] A ; TCALC ; KSCALC ; NDFN ; FN 

C 2 ] A VOLTAGE READING W/O BLACKBODY SOURCE (BBS) = X AND 

[ 3 ] fl VOLTAGE READING W/ BLACKBODY = Y 

[4] fl THIS PROGRAM CALCULATES THE SPECTRAL TRANSMITTANCE 
C 5 ] fl (TAU) AND SPECTRAL EMITTANCE OF A FLAME USING THE 

[6] fl 3.94 MICRON INFRARED SPECTROMETER FILTER AND THE 

[7] fl NEUTRAL DENSITY FILTER ND=3 

[8] fl THIS PROGRAM WILL CALCULATE A CALIBRATION CONSTANT 

[9] fl KS, BASED ON THE VOLTAGE READING FOR THE BLACK BODY 
CIO] fl SOURCE AT 1000 DEGREES C. 

Cll] fl FOUR READINGS OF VOLTAGE ARE REQUIRED: 

C12] fl VC 1 ] =V0 , V[ 2 ] =VBB, V[ 3 ] =VF1 , VC 4 ] =VF2 
C 13 ] NTBB<-0. 051064 fl NT FOR BBS AT 1000 DEGREES C. 

C14] NREF<-4.442E“5 fl NREF FOR C02 FILTER. 

C 15 ] AT<-33 fl CONSTANT ATTENTUATION VALUE. 

C 16 ] SAT<-4 4.6 fl SENSITIVITY FOR ATTENUATION VALUE OF 33. 
C17] 

[18] fl CALCULATION OF THE CALIBRATION CONSTANT KS FROM THE 

[19] fl VOLTAGE READING FOR THE BLACK BODY SOURCE AT 1000 C. 

[20] KSCALC*- ( SATxVBB ) * ( NTBB-NREF ) 

[21] 'KSCALC = ' , (*KSCALC) 

[22] VOLT*-0 , VBB , X , Y 

[23] 

[24] fl CALCULATE THE SPECTRAL TRANSMITTANCE 

[ 25 ] TAU<-( (VOLTC 4 ] xSAT) - CVOLTC 3 ] xSAT) )* ( ( VOLTC 2 ] xSAT) - 

[26] (VOLTC 1] xSAT) ) 

[27] EMM*- 1 -TAU fl CALCULATE THE EMISSIVITY 

[28] fl .964 IS A TRANSMISSION CORRECTION FACTOR, SINCE FLAME 

[29] fl IS <60 ' FROM RAD 

[30] NTF<-NREF+ ( ( SATxVOLT [ 3 ] ) + ( KSCALCxEMMx 0.964)) 

[31] ALN«-®NTF 

[32] 'NTF = ' , (*NTF) 

[33] 

[34] fl USING A CURVEFIT TO THE NT VS TEMPERATURE TABLE, 

[35] fl CALCULATE THE FLAME TEMPERATURE: 

[36] TCALC<-6379+ ( 5072. 53xALN) + ( 2267. 85x(ALN* 2) ) + (659. 18x 

[37] (ALN*3) )+(127.59x(ALN*4) )+(16.23 x(ALN*5))+(1.295x 

[38] (ALN*6) )+(0.0584x(ALN*7) )+(1.129E"3x(ALN*8) ) 

[39] TCALC«-TCALC+273 fl CONVERT TO KELVIN. 

[40] 'A***************************************************' 

[41] 

[42] fl PRINT OUT CALCULATED FLAME TEMPERATURE: 

[43] 'FLAME TEMPERATURE (DEG. K) =',(*TCALC) 

V 
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